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ABSTRACT 

Using the two-fluid formalism with the polytropic approximation, we examine the ax- 
isymmetric stability criterion for a composite system of gravitationally coupled stellar 
and gaseous singular isothermal discs (SIDs). Both SIDs are taken to be infinitely thin 
and are in a self-consistent steady background rotational equilibrium with power-law 
surface mass densities (oc r _1 ) and flat rotation curves. Recently, Lou & Shen (2003) 
derived exact solutions for both axisymmetric and nonaxisymmetric stationary per- 
turbations in such a composite SID system and proposed the D s — criterion of stability 
for axisymmetric perturbations. Here for axisymmetric perturbations, we derive and 
analyze the time-dependent WKBJ dispersion relation to study stability properties. 
By introducing a dimensionless stellar SID rotation parameter D Sl defined as the ra- 
tio of the constant stellar rotation speed V s to the constant stellar velocity dispersion 
a s , one can readily determine the axisymmetric stability D s — criterion numerically 
by identifying a stable range of D s . Those systems which rotate too slow (collapse) 
or too fast (ring fragmentation) are unstable. We found that the stable range of D 2 
depends on the mass ratio 6 of the gaseous SID to the stellar SID and on the square 
of the ratio /3 of the stellar velocity dispersion (which mimics the sound speed) to the 
gaseous isothermal sound speed. Increment of either 6 or (3 or both will diminish the 
stable range of D 2 a . The WKBJ results of instabilities provide physical explanations 
for the stationary configurations derived by Lou & Shen. It is feasible to introduce 
an effective Q parameter for a composite SID system. The closely relevant theoretical 
studies of Elmegreen (1995), Jog (1996) and Shu et al. (2000) are discussed. A study 
of composite partial SID system reveals that an axisymmetric dark matter halo will 
promote stability of composite SID system against axisymmetric disturbances. Poten- 
tial applications to disc galaxies, circumnucleus discs around nuclei of galaxies, and 
protostellar discs are briefly discussed. 

Key words: stars: formation — ISM: general — galaxies: kinematics and dynamics — 
galaxies: spiral — galaxies: structure. 



1 INTRODUCTION 

In the context of galactic disc dynamics, there have been nu- 
merous investigations on criteria of axisymmetric and non- 
axisymmetric disc instabilities (e.g., Binney & Tremaine 
1987). For a single disc of either gaseous or stellar content, 
Safronov (1960) and Toomre (1964) originally derived a di- 
mensionless Q parameter to determine the local stability 
(i.e., Q > 1) against axisymmetric ring-like disturbances. 
Besides a massive dark matter halo, a more realistic disc 
galaxy involves both gas and stars. It is thus sensible to 
consider a composite system of one gas disc and one stel- 



lar disc under the gravitational control of the dark matter 
halo. Theoretical studies on this type of two-component disc 
systems were extensive in the past (Lin & Shu 1966, 1968; 
Kato 1972; Jog & Solomon 1984a, b; Bertin & Romeo 1988; 
Romeo 1992; Elmegreen 1995; Jog 1996; Lou & Fan 1998b, 
2000a, b). In particular, it has been attempted to introduce 
proper definitions of an effective Q parameter relevant to a 
composite disc system for the criterion of local axisymmetric 
instabilities (Elmegreen 1995; Jog 1996; Lou & Fan 1998b). 
The results of such a stability analysis may provide a basis 
for understanding composite disc dynamics (e.g. Lou & Fan 
2000a, b; Lou & Shen 2003) and for estimating global star 
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formation rate in a disc galaxy (e.g., Kennicutt 1989; Silk 
1997). 

There are several reasons that lead us to once again 
look into this axisymmetric stability problem. A few years 
ago, Shu, Laughlin, Lizano & Galli (2000; see also Galli 
et al. 2001) studied stationary coplanar perturbation struc- 
tures in an isopedically magnetized singular isothermal disc 
(SID) without using the usual WKBJ approximation. They 
found exact analytical solutions for stationary (i.e., zero 
pattern speed) configurations of axisymmetric and non- 
axisymmetric logarithmic spiral perturbations. According 
to their analysis, for axisymmetric perturbations with ra- 
dial propagation, a SID with sufficiently slow rotation speed 
will collapse for sufficiently large radial perturbation scales, 
while a SID with sufficiently fast rotation speed will be 
unstable to ring fragmentations for relatively small radial 
perturbation scales (see Fig. 2 of Shu et al. 2000). For ax- 
isymmetric perturbations, the stable regime of SID rotation 
speed may be characterized by a dimensionless rotation pa- 
rameter D, which is the ratio of the constant SID rotation 
speed V to the isothermal sound speed a. The critical values 
of the lowest and the highest D for axisymmetric stability 
can be derived from their marginal stability curves. Shu et 
al. (2000) supported their interpretations by invoking the 
familiar Toomre Q parameter, implying that a local WKBJ 
analysis may still have some relevance or validity to global 
SID perturbation solutions. In their proposed scheme, Shu 
et al. noted that the two critical values of D are fairly close 
to Q = 1 for neutral stability. As they did not invoke the 
WKBJ approximation in deriving the analytic perturbation 
solutions, the Q— criterion and the D— criterion correspond 
to each other well with large wavenumber (ring fragmen- 
tation), while the rough correspondence of the Q— criterion 
and the D— criterion for small wavenumbers (collapse) is un- 
expected where the WKBJ requirement is poorly met. 

Regarding astrophysical applications, these basic prop- 
erties of D— criterion should be properly applied in contexts 
even at qualitative or conceptual levels. For example, for 
majority of currently observed disc galaxies, the typical ro- 
tation parameter D may be sufficiently large to avoid the 
collapse regime, that is, these disc galaxies are rotationally 
supported. Besides other instabilities, the ring fragmenta- 
tion instability is conceptually pertinent to disc galaxies 
such as estimates of global star formation rates (e.g., Jog 
& Solomon 1984a; Kennicutt 1989; Silk 1997; Lou & Fan 
1998b, 2002a, b). Having said these, we should add that the 
collapse regime of D— criterion might be relevant in early 
stages of galactic disc evolution. That is, slow and fast disc 
rotations of proto-galaxies may be responsible for eventual 
bifurcations into different morphologies of galactic systems. 
Furthermore, for central discs (with radii less than a few 
kiloparsecs) surrounding galactic nuclei with supermassive 
black holes (e.g., Lynden-Bell 1969) or for protostellar discs 
around central collapsed cores (e.g., Shu 1977), this D pa- 
rameter may be sufficiently small to induce disc collapses. 
Therefore, both the collapse and ring-fragmentation regimes 



can be of considerable interest in various astrophysical ap- 
plications. 

Recently, we investigated coplanar perturbation config- 
urations in a composite system of two gravitationally cou- 
pled stellar and gaseous SIDs, both taken to be razor thin 
(Lou & Shen 2003). In the ideal two-fluid formalism, we de- 
rived analytical solutions for stationary axisymmetric and 
logarithmic spiral configurations in a composite SID system 
in terms of a dimensionless rotation parameter D s for the 
stellar SID*. By the analogy of a single SID (Shu et al. 2000), 
these stationary configurations should form marginal stabil- 
ity curves (see Fig. 2 and more extensive results of Lou & 
Shen 2003) that separate the stable region from unstable ' 
regions for axisymmetric disturbances. Compared with the 
single SID case, the stable range of D s is reduced owing to 
the presence of an additional gaseous SID and can vary for 
different parameters. 

In order to make convincing arguments for the above 
analogy and our physical interpretations, we perform here 
a time-dependent local stability analysis using the WKBJ 
approximation for the composite SID system (Lou & Shen 
2003). In a proper parameter regime, we demonstrate the 
correctness of interpretations by Shu et al. (2000), and in 
general, we clearly show the validity of using D s parameter 
to demark stable and unstable regimes. To place our analy- 
sis in relevant contexts, we also discuss how the two effective 
Q parameters of Elmegreen (1995) and of Jog (1996) are re- 
lated to our D 3 parameter when other pertinent parameters 
are specified. 

In Section 2, we derive the time-dependent WKBJ dis- 
persion relation for the composite SID system and define 
several useful dimensionless parameters. In Section 3, we 
present the criterion for the SID system being stable against 
all axisymmetric disturbances. Finally, our results are dis- 
cussed and summarized in Section 4. 



2 TWO-FLUID COMPOSITE SID SYSTEM 

Here, we go through the basic fluid equations and derive the 
local dispersion relation in the WKBJ approximation (Lin & 
Shu 1964). For physical variables, we use either superscript 
or subscript s for the stellar SID and g for the gaseous SID. 
As both SIDs are assumed to have fiat rotation curves with 
different speeds V s and V g , we express the mean angular 
rotation speeds of the two SIDs in the forms of 

fi s = Vs/r = a s D s /r, VL g = V g /r = a g D g /r , (1) 

where a s and a g are the velocity dispersion of the stellar 
SID and the isothermal sound speed in the gaseous SID, 
respectively. Dimensionless rotation parameters D s and D g 
are defined as the ratios of stellar SID rotation speed to 
velocity dispersion and gaseous SID rotation speed to sound 



* Subscript s indicates the relevance to the stellar SID. 
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(2) 



speed, respectively. The two epicyclic frequencies are 

Ks ee {{2n a /r)[d{r 2 n B )/dr]} 1/2 = V2Q S , 

K g = {(2n g /r)[d(r 2 n g )/dr]} 1/2 = s/2U a . 

Note that k s and n g are different in general. 

By the polytropic assumption for the surface mass 
density and the two-dimensional pressure (e.g., Binney & 
Tremaine 1987), the fluid equations for the stellar disc, in 
the cylindrical coordinates (r, ip, z) at the z = plane, are 

^ s + I^ (rS v ) + lA (sY) = 0i (3 ) 



dt 



r dr 



du s s du s f du 3 
dt dr r 2 dip 



1 d t 2^ d(f> , . 
E^^ KE) -^' (4) 



dj 



, dj s , f dj 



1 d 



dt 



dr 



dp 



E s dp' 



dp 



(5) 



where E s is the stellar surface mass density, u s is the radial 
component of the bulk fluid velocity, j" is the z— component 
of the specific angular momentum, a a is the stellar velocity 
dispersion (or an effective "isothermal sound speed"), a 2 T, s 
stands for an effective two-dimensional pressure in the poly- 
tropic approximation, and (f> is the total gravitational poten- 
tial perturbation. For the corresponding fluid equations in 
the gaseous disc, we simply replace the relevant superscript 
or subscript s with g. The two sets of fluid equations are 
coupled by the total gravitational potential 4> through the 
Poisson integral 

-GT,(r',i:,t)r'dr' 



(r, <p,t)= f dip 



o [r' 2 + r 2 — 2rr' cos(ip — <£>)] 



(6) 



where E = E s +E 9 is the total surface mass density. Here the 
stellar and gaseous SIDs interact mainly through the mutual 
gravity on large scales (Jog & Solomon 1984a, b; Bertin & 
Romeo 1988; Romeo 1992; Elmegreen 1995; Jog 1996; Lou 
& Fan 1998b, 2000a,b; Lou & Shen 2003). 

Using the above equations for the steady background 
rotational equilibrium (indicated by an explicit subscript 0) 



with uq 



0, Q s = jo/r , and fi 9 = Jo/ r > we obtain 



Eg = a 2 (l + D 2 )/[27TGr(l + 5)] , 

X 9 =a 2 g (l + D 2 )S/[2nGr(l + 5)] , 

and an important relation^ 

a 2 (D 2 3 + l)=a 2 g (D 2 g + l) , 

where S = Eq/Eq is the SID surface mass density ratio, 
P ee a 2 /a 2 stands for the square of the ratio of the stellar 



(7) 



(8) 



T This implies a close relation among the two SID rotation speeds, 
the stellar velocity dispersion and the gas sound speed. In this as- 
pect, it is different from the usual local prescription of a composite 
disc system (e.g., Jog &; Solomon 1984a, b; Bertin & Romeo 1988; 
Elmegreen 1995; Jog 1996; Lou & Fan 1998b). Once the ratio of 
a 2 to a 2 is known, the rotation speeds of the composite SID sys- 
tem can be determined through that of the stellar disc, i.e. D 2 , as 
the two rotation speeds are coupled dynamically. In the limiting 
regime of a <C V and j3 ^> 1, V s and V g are comparable with V g 
slightly larger by an amount of ~ (a 2 — a 2 )/(2V a ). In this sense, 



velocity dispersion to the gas sound speed. In our analysis, 
we specify values S and ft to characterize different composite 
SID systems. For late-type spiral galaxies, gas materials are 
less than the stellar mass with 5 < 1. For young proto disc 
galaxies, gas materials may exceed stellar mass in general 
with 5 > 1. Thus, in our analysis and computations, both 
cases of S < 1 and S > 1 are considered. As the stellar 
velocity dispersion is usually larger than the gas sound speed 
(a 2 > a 2 ), we then have /3 > 1. When /3 = 1, we have 
D s = D g by condition @, which means the two SIDs may 
be treated as a single SID (Lou & Shen 2003). 

For small perturbations denoted by subscripts 1, we ob- 
tain the following linearized equations 



<9E? 1 d . _ s , n 9Ef Eg dj{ 

~5T + --5-( rE o^i) + 1" 2 H 

dt r dr dip r z dip 







dt s dip 



ji d ( 2 E? 



0) 



(10) 



dj! , k 2 s dj{ d ( 2 Ef \ 

for the stellar SID, and the corresponding equations for the 
gaseous SID with subscript or superscript s replaced by g 
in equations (9) — (11). The total gravitational potential per- 
turbation is 



t/)i(r,ip,t) = fdip 



-G(Ef + T,f)r'dr' 



(12) 



the assumption k s 



is reasonable. 



/o [r' 2 + r 2 — 2rr' cos(ip — ip)] 1 ^ 2 

In the usual WKBJ approximation, we write the coplanar 
axisymmetric perturbations with the periodic dependence of 
exp(ikr + icut) where k is the radial wavenumber and to is the 
angular frequency. With this dependence in the perturbation 
equations, we obtain the local WKBJ dispersion relation for 
the composite SID system in the form of 

(lo 2 — k 2 — k 2 a 2 + 27rGjfcjEQ) 

x (o; 2 -Kg-fc 2 aj + 27rG|fc|Eg) (13) 

- (27rG|fc|Eo)(27rG|fc|Eg) = . 

Jog & Solomon (1984a) have previously derived the similar 
dispersion relation with V s = V g and thus k s = n g . The 
stability analyses of Elmegreen (1995) and Jog (1996) start 
with the same dispersion relation of Jog & Solomon (1984a). 
The WKBJ approach of Jog & Solomon (1984a) is generally 
applicable to any locally prescribed properties of a two- fluid 
disc system. In earlier theoretical studies including that of 
Jog & Solomon (1984a), it is usually taken that the rota- 
tion speeds of the two fluid discs are equal. Our model of 
two coupled SIDs here is self consistent globally with a cen- 
tral singularity and with different SID rotation speeds in 
general. When one applies the local WKBJ analysis to our 
composite SID system, the procedure is identical to that of 
Jog & Solomon (1984a) but with different local SID speeds 
V s and V g and hence different epicyclic frequencies k s and 
K g . Therefore, our dispersion (13) appears strikingly similar 
to that of Jog & Solomon (1984a) yet with different k s and 
Kg. In short, if one prescribes two different disc speeds and 
thus two different epicyclic frequencies in the analysis of Jog 
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& Solomon (1984a), the resulting dispersion relation should 
be our equation (13). 

In our analysis, the rotation speeds of the stellar and 
the gaseous SIDs differ in general. In particular, since a 2 
is usually higher than a 2 , it follows from equations (1) and 
(8) that Q.s < £l g . That is, the stellar disc rotates some- 
what slower than the gaseous disc does. Such phenomena 
have been observed in the solar neighborhood as well as 
in external galaxies (K. C. Freeman, 2003, private commu- 
nications). This is related to the so-called asymmetric drift 
phenomena (e.g., Mihalas & Binney 1981) and can be under- 
stood in terms of the Jeans equations (Jeans 1919) originally 
derived by Maxwell from the collisionless Boltzmann equa- 
tion. In essence, random stellar velocity dispersions produce 
a pressure-like effect such that mean circular motions be- 
come slower (Binney & Tremaine 1987). In our treatment, 
this effect of stellar velocity dispersions is modelled in the 
polytropic approximation (see eqns (2) and (8)). Biermann 
& Davis (1960) discussed the possibility that the mean ro- 
tation speeds of the stars and gas in the Galaxy may be dif- 
ferent and reached the opposite conclusion that a gas disc 
should rotate slower than a stellar disc does. However, their 
analysis based on a virial theorem that excludes the effect of 
stellar velocity dispersions. In other words, the stress tensor 
due to stellar velocity dispersions in the Jeans equation has 
been ignored in their considerations. 



3 AXISYMMETRIC STABILITY ANALYSIS 

We study the axisymmetric stability problem based on 
the WKBJ dispersion relation. Similar to previous analy- 
sis (Toomre 1964; Elmegreen 1995; Jog 1996; Lou & Fan 
1998b), we also derive an effective Q parameter. This Q e ff 
is expected to depend on D 2 once 8 and /3 are known. 
For the convenience of analysis, we define notations 

H\ = k 2 + k 2 a 2 — 27rG|fc|Eo , 



(14) 



H 2 = Kg + k 2 a 2 g — 27tG|A:|Eq , 

Gi = 2 7 rG|/c|E i , 

G 2 = 27rG|fc|Eg , 

to express dispersion relation l|13|l as an explicit quadratic 
equation in terms of u> 2 , namely 

uj 4 - {Hi + H 2 )uo 2 + (HiH 2 - G1G2) = . 

The two roots* w\ and uPl of equation 115H are 

<A(k) = \{(Hi + H 2 ) 

± {(Hi + H 2 ) 2 - 4(HiH 2 - GiG 2 )] 1/2 } 
The uj\ root is always positive, which can be readily proven. § 



(15) 



(16) 



t The determinant A = (Hi — H 2 ) 2 + 4GiG 2 > 0. 

§ This conclusion holds true if Hi + H 2 > 0. If Hi + H 2 < 0, 

then HiH 2 — GiG 2 < and the conclusion still holds. 



We naturally focus on the (J 2 , root, namely 
w-(*) = i{(ffi + ffa) 



(17) 



- [(Hi + H 2 y - A(HiH 2 - GiG 2 )] ' } , 

to search for stable conditions that make the minimum of 
Lui(k) > for all k. 

By expanding the to'L(k) root in terms of definition 1141 . 
we obtain an expression involving parameters n 2 , k, 2 , a 2 , a 2 , 
Eq, Eq, and k. For a composite SID system with power-law 
surface mass densities, flat rotation curves, and the equilib- 
rium properties of k s , re 9 , Eg and Eq , equation 117H can also 
be expressed in terms of four parameters a 2 , D 2 , 5 and /?, 
namely 

u-(K) = -^{A 2 K 2 + AiK + A - p 1/2 ) , (18) 

where K = \k\r is the dimensionless local radial wavenumber 
and the relevant coefficients are defined by 

A 2 = 1 + 1/p , 



Ai = -(1 + D 2 ) , 
A = 4D 2 + 2{l-l//3) 



(19) 



with 



B 4 K A + B3K 3 + B 2 K 2 + BiK + Bo , (20) 



B A = (1 - 1//3) 2 , 

B 3 = 2(1 + D 2 )(l - l/0)(5 - !)/(! + S) , 



B 2 ee (D 2 - 1 + 2/(5)(D 2 + 3 - 2/p) , 
Bi = 4(1 + D 2 )(l - 1/(5)(1 - 8)/(l + 8) , 
Bo = 4(1 - 1/P) 2 . 



(21) 



As a 3 has been taken out as a common factor of equation 
H18L the most relevant part is simply 

(22) 



A 2 K 2 +AiK + A -p 1/2 



that involves the stellar rotation parameter D 2 and the two 

ratios S and ft. We need to determine the value of K m ^ n at 

which ^(K) takes the minimum value, in terms of D 2 . We 

then substitute this J^ mm into expression 12211 and derive 

the condition for D 2 that makes this minimal uj 2 - > 0. 

3 mm — 



3.1 The D s — criterion in the WKBJ regime 

To demonstrate the stability properties unambiguously, and 
to confirm our previous interpretations for the marginal sta- 
bility curves in a composite SID system (Lou & Shen 2003), 
we first present some numerical results that contain the same 
physics and are simple enough to be understood. 

According to equation 1181 . uj'L is a function of rota- 
tion parameter D 2 and wavenumber K. We then use K as 
the horizontal axis and D 2 as the vertical axis to plot con- 
tours of uit numerically. A specific example of 8 = 0.2 and 
(3 = 10 is shown in Fig. 1. The shaded regions are where (JL 
takes on negative values and the blank region is where (JL 
takes on positive values. In this manner, we clearly obtain 
marginal stability curves (solid lines) along which lo 2 . = 0. 
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Figure 1. Contours of as a function of if and D 2 with <5 = 0.2 
and /3 = 10. The uP_ > region is blank and the < regions 
are shaded. 

Physically, the shaded regions of negative w 2 . values are un- 
stable, while the blank region of positive uj 2 _ values is stable. 
From such contour plots of lu 2 _, one can readily determine 
the specific range of stellar rotation parameter D 2 such that 
the composite SID system is stable against axisymmetric 
disturbances. We therefore confirm our previous interpreta- 
tions for the marginal stability curves for a composite SID 
system (Lou & Shen 2003). For a comparison, Fig. 2 shows 
the corresponding marginal stability curves with 8 = 0.2 
and P = 10 that we derived previously (Lou & Shen 2003) 
as the stationary axisymmetric perturbation configuration 
where a is a dimensionless effective radial wavenumber (see 
also Shu et al. 2000) . The apparent difference in Figs. 1 and 
2, mainly in the collapse regions, arises because the results 
of Fig. 2 are exact perturbation solutions without using the 
WKBJ approximation. 

Now in the WKBJ approximation, the stable range of 
D 2 can be read off from Fig. 1. For the case of 8 = 0.2 
and P — 10, the composite SID system is stable from 
D 2 = 0.1193 at K = 0.5121 to D 2 = 4.1500 at K = 3.9237. 
In comparison with the single SID case (see Fig. 2 of Shu 
ct al. 2000), this stable range of D 2 diminishes. With other 
sets of 8 and /3 parameters, we can derive similar diagrams 
of Fig. 2 as have been done in our recent work, for examples, 
Figs. 7 - 10 of Lou & Shen (2003). In the WKBJ approxi- 
mation, we show variations with 8 and j3 in Figs. 3 and 4. 
When 8 is larger, the system tends more likely to be unsta- 
ble for small (3. In reference to Fig. 10 of Lou & Shen (2003) 
with 8 — 1 and /3 = 30, we plot contours of u>_ in Fig. 5 
in the WKBJ approximation. In this case, the stable range 
of D 2 disappears, that is, the composite SID system can- 
not be stable against axisymmetric disturbances. We note 
again that the collapse region with small wavenumber devi- 
ates from the exact result of Fig. 10 of Lou & Shen (2003) , 
due to the obvious limitation of the WKBJ approximation. 
For the unstable region of large wavenumber, which we la- 
belled as ring fragmentation previously (Lou & Shen 2003) , 
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stable oscillations 
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Figure 2. The marginal stability curves of D 2 versus the dimen- 
sionless effective radial wavenumber a for m = 0, <5 = 0.2, and 
(i = 10. While <5 remains small, a larger ft lowers the ring fragmen- 
tation boundary. It is then easier for the SID system to become 
unstable in the form of ring fragmentation but with a larger a. 

Fig. 5 here and Fig. 10 of Lou & Shen (2003) show good 
correspondence as expected. For the purpose of comparison, 
Table 1 lists several sets of 8 and j3 and the corresponding 
stable range of D 2 , selected from our stationary perturba- 
tion configurations studied earlier (see Figs. 5-10 of Lou & 
Shen 2003), using both the WKBJ approximation here and 
the marginal stability curves of Lou & Shen (2003) . 

It is not surprising that there exist one lower limit and 
one upper limit to bracket the stable range of D 2 . We now 
provide interpretations for the single SID case. In the usual 
WKBJ analyses, the surface density E is prescribed and 
hence independent of the epicyclic frequency k. In our global 
SID formalism, the surface mass density and the epicyclic 
frequency are coupled through the rotation parameter D. 
Both E and k decrease with decreasing D but with different 
rates. Specifically, the surface mass density E is proportional 
to 1 + D 2 (see eq. (7)) while k is proportional to D (see eq. 
(2)). For either large or small enough D, the ratio of k to 
E may become less than the critical value (sound speed is 
constant by the isothermal assumption), that is, the Q pa- 
rameter becomes less than unity for instabilities. 

We stress that ring fragmentation instabilities occurring 
at large wavenumbers are familiar in the usual local WKBJ 
analysis (Safronov 1960; Toomre 1964). While collapse in- 
stabilities at small wavenumbers are known recently for the 
SID system (Shu et al. 2000; Lou 2002; Lou & Shen 2003). 
The physical nature of the new "collapse" regime is funda- 
mentally due to the Jeans instability when the radial pertur- 
bation scale becomes sufficiently large. 1 By the conservation 

» One might wonder how a large-scale disc Jeans instability can 
be revealed in a local WKBJ analysis although somewhat crudely. 
Perhaps, the best analogy is the WKBJ dispersion relation for 
spiral density waves in a thin rotating disc (Lin & Shu 1966, 
1968). This WKBJ dispersion relation is quadratic in the radial 
wavenumber |fe| and can be solved for a larger |fe| and a smaller 
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Figure 3. A contour plot of ui 2 _ as a function of K and D 2 with 
5 = 1 and fi = 5. The lo 2 _ > region is blank and the ui 2 ^ < 
region is shaded. 
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Figure 5. A contour plot of uP_ as a function of K and D 2 with 
8 = 1 and j3 = 30. The u?_ > region is blank and the io 2 _ < 
region is shaded. 




co 2 >0 

a' 1 1 1 1 1 1 1 

012345678 

K 

Figure 4. A contour plot of uP_ as a function of if and D 2 with 
5 = 5 and /3 = 3. The u>£. > region is blank and the oj^ < 
region is shaded. 

of angular momentum, SID rotation plays a stabilizing role 
to modify the onset of Jeans collapse (e.g., Chandrasekhar 
1961). Hence, there appears a critical D 2 below which exists 
the collapse regime. Of course, too large a D 2 causes the 
other type of ring fragmentation instabilities to set in. For 
most disc galaxies currently observed, the rotation parame- 
ter D 2 is typically large enough to avoid such collapses. 

Elsewhere, the collapse regime studied here may have 



\k\ corresponding to short- and long- wave branches respectively 
(e.g., Binney & Tremaine 1987). The former is naturally consis- 
tent with the WKBJ approximation, while the latter, inconsistent 
with the WKBJ approximation, turns out to be a necessary key 
ingredient in the swing process (Goldrcich & Lyndcn-Bcll 1965; 
Toomre 1981; Fan & Lou 1997) as revealed by analytical and nu- 
merical analyses. In the end, one needs to perform global analyses 
to justify the WKBJ hint for instabilities of relatively large scales. 



relevant astrophysical applications, e.g., in circumnuclear 
discs around nuclei of galaxies and protostellar discs around 
protostellar cores etc., where D 2 may be small enough (ei- 
ther a high velocity dispersion a s or a low rotation speed 
V s ) to initiate collapses during certain phases of disc system 
evolution. In fact, for a low rotation speed V s , we note from 
the perspective of evolution, that those proto disc galaxies of 
low rotation speeds in earlier epochs, will Jeans collapse and 
lead to other morphologies, while those proto disc galaxies 
of sufficiently fast rotation speeds in earlier epochs will grad- 
ually evolve into disc galaxies we observe today. Meanwhile, 
ring fragmentation instabilities are thought to be relevant to 
global star formation in a disc galaxy - another important 
aspect of galactic evolution. 

The addition of a gaseous SID to a stellar SID will 
decrease the overall stability against axisymmetric distur- 
bances at any wavelengths in general, that is, the stable 
range of D 2 shrinks. However, for the two types of instabil- 
ities, namely, the collapse instability at small wavenumbers 
and the ring-fragmentation instability at relatively large 
wavenumbers, the gravitational coupling of the two SIDs 
play different roles. For instance, a composite SID system 
tends more likely to fall into ring fragmentations but less 
likely to collapse. For either larger 8 or larger f3, this ten- 
dency becomes more apparent. 

This trend can be understood physically. For the ring- 
fragmentation regime, suppose the stellar SID by itself alone 
is initially stable and a gaseous SID component is added to 
the stellar SID system in a dynamically consistent manner. 
The gaseous SID may be either stable or unstable by itself 
alone. By itself, the gaseous SID tends to be more unstable 
if the sound speed a g becomes smaller or £q becomes larger 
according to the definition of the Q parameter (valid for 
the ring fragmentation instability) , which means larger j3 or 
larger 5, respectively. Thus, with either (5 or 8 or both being 
larger, a composite SID system tends to be more vulnerable 
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10 


0.1193(0.7602) 
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30 
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1 


1.5 


0.0596(0.6556) 
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0.0039(0.1418) 


1.7397(1.6718) 
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5 


O.OOll(-) 


0.5696(0.5317) 



Table 1. The overall stable ranges of D 2 against axisymmetric 
disturbances at any wavelengths for different sets of 8 and fi. The 
values in parentheses are derived by Lou & Shen (2003) without 
the WKBJ approximation and are preferred to describe the col- 
lapse region. The overall stable range of appears to decrease 
when cither <5 or /3 increases. The stable range of D 2 S may disap- 
pear for a case like 8 = 1 and f3 = 30. When (3 = 1, the problem is 
independent of 8, that is, <5 can be arbitrary (Lou & Shen 2003). 



to ring fragmentations. The interpretation for the collapse 
regime is associated with the dynamical coupling of surface 
mass density E and epicyclic frequency k that leads to the 
following conclusion, namely, a lower velocity dispersion of 
the gas component seems to prevent collapse. This should 
be understood from condition (8): for larger f3, the rotation 
speed V g will exceeds V s by a larger margin and helps to 
prevent an overall SID collapse. 

As already noted, the WKBJ approximation becomes 
worse in the quantitative sense when dealing with small 
wavenumbers where the collapse regime exists. But the sta- 
bility properties of a SID can be qualitatively understood 
through the WKBJ analysis. The exact perturbation analy- 
sis of our recent work (Lou & Shen 2003) has shown that the 
collapse regime diminishes with increasing values of 5 and (3. 
Therefore, once the interpretations are justified physically, 
the exact procedure of Lou & Shen (2003) should be adopted 
to identify the relevant stable range of D 2 of a composite 
SID. As much as Shu et al. (2000) have done for a single SID 
analytically, Lou & Shen (2003) did the same for a composite 
SID system with exact solutions to the Poisson integral. The 
stationary axisymmetric perturbation configurations were 
derived to determine the marginal stability curves. These 
marginal stability curves have been proven to have the same 
physical interpretation as the WKBJ marginal curves ob- 
tained here. 



3.2 Effective Q parameters 

It is natural to extend the concept of the Q parameter for 
a single disc to an effective Q parameter for a composite 
disc system for the purpose of understanding axisymmetric 
stability properties. In the following, we discuss effective Q 



parameters for a composite SID system in reference to the 
works of Elmegreen (1995) and Jog (1996). There are sev- 
eral points worth noting in our formulation. First, n s and 
K>g are related to each other but are allowed to be different 
in general. Secondly, the background surface mass densities 
are related to SID rotation speeds through the polytropic 
approximation and the steady background rotational equi- 
librium. For the convenience of notations, the effective Q pa- 
rameters introduced by Elmegreen (1995) and by Jog (1996) 
are denoted by Qe and Qj, respectively. 



3. 2. 1 The Q e parameter of Elmegreen 

Starting from equation (18) or (22), we define an effective 
Q e parameter in a composite SID system following the pro- 
cedure of Elmegreen (1995) but with two different epicyclic 
frequencies k s and n g . The minimum of u) 2 _ is given by 

„2 

+ A^K min + A - p 1/2 ] 

(23) 



2r 



2 \AiK, 
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mm 
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A'lK-min — AlK„ 



A 



0(Ql 



1) 



x / 2 - A 2 Kl 



(24) 



A\K m i n 

and p takes on the value at K = K min . As K min depends 
on D 2 , 8 and /3, so does Q%. For specified parameters 8 
and P of a composite SID system, the parameter Q%, giv- 
ing the criterion of axisymmetric stability (i.e. stable when 
Q 2 E > 1), corresponds to the rotation parameter D 2 of the 
stellar disc. One purpose of the present analysis is to estab- 
lish the correspondence between the axisymmetric stability 
condition and the marginal stability curves of axisymmet- 
ric stationary perturbation configuration derived recently by 
Lou & Shen (2003) for a composite SID system. The per- 
tinent results for a single SID can be found in Shu et al. 
(2000). 

Our task now is to find the solution K m in- As k s / 
K, g in general, our computation turns out to be somewhat 
more complicated or involved than that of Elmegreen (1995) 
who took k, s = Kg. It becomes more difficult to solve this 
algebraic problem of order higher than the third analytically. 
However, if one multiplies both sides of equation (17) by the 
positive root lo + , a simple expression appears 

G1G2 ■ (25) 

Instead of minimizing uj 2 _, we solve for K c in terms of D 2 
at which W = uj+uj 2 - attains the minimum value. This K c 
must satisfy the following equation 

dK dK ' 

By substituting Hi, H2, Gi and Gi into equation (26) and 
regarding 8 and f3 as specified parameters, equation (26) can 
be reduced to a cubic equation of K involving D 2 , namely, 

K 3 + aK 2 + bK + c = , (27) 

where 



W = w+w- =H 1 H 2 



a = -3(D 2 + l)(l + p5)/[4(l + 8)} 
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b = D 2 s (f3 + i)+{3~l , 



c = -(D 2 + 1)[(3D 2 (1 + 8) +(3- 1]/[2(1 + 8)] . 

In most cases " , equation (27) has only one real solution in 
the form of 

K c ={x- q/2) 1 ' 3 + (-x - g/2) 1/3 - a/3 , (28) 

where x = (g 2 /4 + p 3 /27) 1/2 , p = b-a 2 /3 and q = 2a 3 /27- 
ab/3 + c. For a single real root, it is required that q 2 /4 + 
p 3 /27 > . Through numerical computations, we find that 
this condition is met for most pairs of 8 and /3. 

We then regard this K c as an estimator** for K m i n to 
determine Q%. As already known in our analysis (Lou & 
Shen 2003), when the ratio f3 is equal to 1, the properties 
of a composite SID system have something in common with 
those of a single SID. Moreover, the stationary axisymmetric 
configuration becomes the same as a single SID. We expect 
that the present problem should reduce to a single SID case 
when = 1. We consider below the special case of j3 = 1. 

With (3 = 1 in equation 1)27^ . K c can be determined in 
terms of D 2 without involving 8 at all. We then use K c as 
K m i n in the definition of Q 2 E , namely, condition 12 H . One 
can readily verify that Q% > 1 for 3 - 2^2 < D 2 < 3 + 2^2; 
Qe < 1 for D 2 > 3 + 2^2 or D 2 < 3 - 2^2; and Q% = 1 
for D 2 = 3 ± 2^2. The range of D 2 that makes Q% > 1 
is 0.1716 < D 2 < 5.8284, which is exactly the same as the 
Q estimator for the single SID of Shu et al. (2000). Based 
on estimates of Q parameter alone, Shu et al. (2000) made 
physical interpretations for their axisymmetric stationary 
perturbation solution curves for a single SID. According to 
their interpretation, the maximum D 2 = 0.9320 of the col- 
lapse branch and the minimum D 2 = 5.410 of the ring- 
fragmentation branch (see Fig. 2 of Shu et al. 2000) encom- 
pass the required range of D 2 for stability against axisym- 
metric disturbances. As the formulation of Shu et al. (2000) 
is exact in contrast to the WKBJ approximation, the sta- 
ble range of D 2 slightly deviates from that determined by 
the effective Qe parameter. And for the upper bound of D 2 , 
with a larger wavenumber (or a smaller wavelength) , the two 
approaches lead to consistent results as expected. 

In our recent analysis (Lou & Shen 2003), we obtained 
curves for stationary axisymmetric perturbations given var- 
ious combinations of 8 and j3 for a composite SID system. 
These curves differ from those for a single SID and the sta- 
ble range of D 2 varies for different combinations of 8 and 
(3. It is found that this stable range of D 2 may be signif- 
icantly reduced for certain values of 8 and f3, as shown in 
Fig. 7 — 9 of Lou & Shen (2003) and in the analysis of Sec- 
tion 3.1. Moreover, when 8 and f3 are sufficiently large, the 

II When k s = Kg, Bertin & Romeo (1988) give the proof that 
as long as > 0.0294 and 8 > 0.172, there exists only one 
minimum of uP_ . When there are more than one local minima for 
uP_ , one should choose K c for the smallest minimum of uP_ . 
** While K c is obtained as the value at which uP_ multiplying lJ^_ 
reaches the minimum. It is valid to use K c to determine whether 
Q E > 1 or Q E < 1 by realizing that u)^_ remains positive and the 
critical condition Q E = 1 is equivalent to oj^oj 2 = 0. 




°D 1 2 3 4 Q 5 | 7 B 

Ds 

Figure 6. Curves of Q E versus D 2 which intersect the horizontal 
line Q E = 1. For each curve, the two points of intersection at 
Q E = 1 give the stable range of D 2 for a composite SID against 
axisymmetric disturbances. The fi = 1 case gives the same result 
of a single SID. 

stable range of D 2 may no longer exist, as shown in Fig. 10 
of Lou & Shen (2003) and Fig. 5 here. For these cases, a 
composite SID system becomes vulnerable to axisymmetric 
instabilities. 

For 8 — 0.2 and (3 — 10, as shown in Fig. 6 of Lou & 
Shen (2003) and also Fig. 2 of this paper, we search for the 
approximate stable range of D 2 with the WKBJ dispersion 
relation 118H for a composite SID system. By using the Q 2 E 
defined by equation 12411 . one can directly plot the curve of 
Q% versus D 2 as shown in Fig. 6, from which one can identify 
the approximate stable range of D 2 from ~ 0.1193 to ~ 
4.1500 where Q 2 E = 1. Several other curves with different 8 
and f3 are also shown in Fig. 6. 

One thing that should be emphasized is that the Q 2 E 
parameter as defined by equation 1241 must be positive, 
otherwise it may occasionally lead to incorrect conclusions 
regarding the sign of uj 2 . as defined by equation (23). More 
detailed investigations indicate that it is only applicable to 
cases where a composite SID system can be stabilized in a 
proper range of D 2 . For special cases that cannot be stabi- 
lized at all, the Q 2 E thus defined may not be relevant. 

3.2.2 The Qj parameter of Jog 

The Q E parameter can be derived by solving cubic equation 
1271 1 for K. It may happen that equation 12711 gives three real 
solutions. In order to obtain an analytical form of K m i n , 
one must identify the absolute minimum of uj 2 _ . This proce- 
dure can be cumbersome. Alternatively, Jog (1996) used a 
seminumerical way to define another effective Q parameter, 
referred to as Qj here, that is different from Q E defined by 
equation 1241 . 

According to solution (17), u'L will be positive or neg- 
ative if H1H2 — G1G2 > or < 0. The critical condition of 
neutral stability is thus H1H2 — G1G2 = 0, which can be 



© 0000 RAS, MNRAS 000, 000-000 



Axisymmetric Stability Analysis 9 



a 



5=0.2 


— P=1 
p=1 .5 


,7/ V C\ 
'// N 


P=10 _ 
-- - p=30 


( 7 v. v\ 

- // % 




I 




I 




1 2 3 4 

D 

s 


5 6 



defined as 



Figure 7. Curves of Q 2 j as function of Di 



~> 2 , with 5 = 0.2 but for 
different j3 values. The stable ranges of are determined by 
the two points of intersection with the horizontal line Q 2 j = 1. 
Stability corresponds to Q 2 j > 1 when D 2 falls between the two 
intersection points. The two intersection points move to the left 
as increases. 
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(see eqns. (5) and (6) of Jog 1996). It follows that Q 2 j < 1 
or > 1 correspond to instability or stability, respectively. 
Therefore, for a given set of 8, /J and D 2 , one can numeri- 
cally determine the value of K m in- By inserting this K m i n 
into equation lloOl , one obtains the value of Q 2 j for any given 
set of {S,P,D 2 }. It is then possible to explore relevant pa- 
rameter regimes of interest. Jog (1996) introduced three pa- 
rameters, namely, the two Toomre Q parameters for each 
disc and the surface mass density ratio of the two discs, to 
determine Q 2 j. Equivalently, we use three dimensionless pa- 
rameters 8, (3 and D 2 instead. Figs. 7 and 8 show variations 
of Q 2 with D 2 . Similar to Q%, the value of Q 2 j > 1 when D 2 
falls in the same ranges. When 8 is fixed and (3 increases, the 
left and right bounds of the stable D 2 range move towards 
left together, while the width of the stable range appears to 
decrease. We have revealed the same trend of variations in 
our recent work (Lou & Shen 2003). It is possible for the 
left bound to disappear when /3 is sufficiently large, and it 
is easier to achieve this when 8 becomes larger. For a suf- 
ficiently large 8, the increase of /3 may completely suppress 
the stable range of D 2 as shown in Fig. 5 for 5=1 and 
(3 = 30. 

In comparison to Qe, the Q.j definition 1301 1 is valid 
for all parameter regimes and avoids improper situations for 
unusual 8 and f3 values. However, to find the value of Q 2 j 
one must perform numerical exploration for each given D 2 , 
while for Q%, the procedure is analytical as long as there 
exists a stable range of D 2 . Regardless, the critical values of 
D 2 found by Q% and Q 2 are equivalent. 



Figure 8. Curves of Q 2 j as function of D 2 with <5 = 0.2 but 
for different j3 values. These curves intersect with the horizontal 
line Q 2 j = 1. For a larger 8 than that of Fig. 7, the two points 
of intersection move much faster as /3 increases. And the stable 
range of D 2 no longer exists when /3 is sufficiently large. 



cast into the form of 

27rGfcEg 



+ 



2nGkY, 9 



= 1 



One can now define a function T such that 
j. _ 2 7 rG'fcSg + 27rGfcSg 



K(D 2 S + 1)/(1 + 8) Kf3(D 2 + 1)8/(1 + 8) 



(29) 



2D 2 + K 2 



+ 



2[/3(D 2 + !)-!] +K 2 



where expressions (2), (7), (8) and K = \k\r are used. We 
then search for K m i„ at which a composite SID system be- 
comes hardest to be stabilized, that is, when ld 2 _ reaches 
the minimum. The effective Q parameter, Qj, can thus be 



3.3 Composite partial SID system 

In disc galaxies, there are overwhelming evidence for the ex- 
istence of massive dark matter halos as inferred from more 
or less flat rotation curves. To mimic the gravitational effect 
of a dark matter halo, we include gravity terms d<&/dr and 
d<&/dtp in the radial and azimuthal momentum equations (4) 
and (5), respectively, where $ is an axisymmetric gravita- 
tional potential due to the dark matter halo. It is convenient 
to introduce a dimensionless parameter F = </>/(</> + $) for 
the fraction of the disc potential relative to the total po- 
tential (e.g., Syer & Tremaine 1996; Shu et al. 2000; Lou 
2002). The background rotational equilibrium is thus mod- 
ified by this additional <E> accordingly. As before, we write 
Jl s = a s D s /r, Q g = a g D g /r, and k s — \/~2Q, s , K g = \/2fi 9 . 
The equilibrium surface mass densities now become 

a 2 (l + D 2 )F " 2n 1 n2 ^ 



So 



yg 

^0 



a- g (l + D 2 )FS 



2-kGt(1 + 8)' u 2nGr(l + S) ' 

where < F < 1. In our perturbation analysis, the dy- 
namic response of this axisymmetric massive dark matter 
halo to coplanar perturbations in a composite SID system is 
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ignored™ (Syer & Tremaine 1996; Shu et al. 2000; Lou 2002). 
We thus derive similar perturbation equations as those ol a 
full SID case (F = 1) and consequently the similar disper- 
sion relation as equation (13) yet with modified equilibrium 
properties. Now following the same procedure of WKB J per- 
turbation analysis described in Section 3.1, we can plot con- 
tours of uj 2 _ as a function of rotation parameter D 2 and 
wavenumber K for a composite partial SID system. The 
case of F = 1 corresponds to a composite full SID system 
that has been studied in Section 3. When F becomes less 
than 1 (i.e., a composite partial SID system), which means 
the fraction of the dark matter halo increases, the stable 
range of D 2 becomes enlarged, as shown in Figs. 9 and 10 
for 8 = 0.2 and /3 = 10. From these contour plots, it is clear 
that the introduction of a dark matter halo tends to stabilize 
a composite partial SID system. For late-type disc galaxies, 
one may take F = 0.1 or smaller. Such composite partial 
SID systems are stable against axisymmetric disturbances 
in a wide range of D 2 . Moreover, those composite full SID 
systems that are unstable may be stabilized by the presence 
of a dark matter halo, as shown by the example of Fig. 11 
for 5 = 1, /3 = 30 and F = 0.5 (this case is unstable for a 
composite full SID system). 



4 DISCUSSIONS AND SUMMARY 

Our model of a composite (full or partial) SID system is an 
idealization of any actual disc systems, such as disc galaxies, 
circumnuclear discs surrounding nuclei of galaxies, or pro- 
tostellar discs. One key feature of the SID model is the flat 
rotation curve. One major deficiency in galactic application 
is the central singularity that sometimes even gives rise to 
conflicting theoretical interpretations (see Shu et al. 2000). 
The isothermality is yet another simplifying approximation. 
In galactic applications, the common wisdom is to introduce 
a bulge around the center to avoid the central singularity. 
By a local WKBJ stability analysis of axisymmetric pertur- 
bations, this model can be utilized as usual with qualifica- 
tions. We now discuss the WKBJ results in the context of 



tt Strictly speaking, gravitational perturbation coupling between 
a thin, less massive, rotating disc and a more massive dark mat- 
ter halo together with proper matching and boundary conditions 
should lead to an increasing number of modes that may or may 
not be stable. We simplify by noting numerical simulations have 
indicated that a typical galactic dark matter halo has a relatively 
high velocity dispersion of the order of several hundred kilometers 
per second. The rationale of the simplification is that perturba- 
tion coupling between those in the SID system and those in the 
dark matter halo becomes weak for a high velocity dispersion in 
the dark matter halo. For example, for a non-rotating spherical 
dark matter halo, effective "acoustic" and "gravity" modes (i.e., 
p— modes and g— modes in the parlance of stellar oscillations; Lou 
1995) can exist within; these p— and g— modes are expected to 
be Landau damped in a collisionless system. However, the stabil- 
ity of "interface modes" between a thin rotating disc and a dark 
matter halo remains to be investigated. 
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Figure 9. Contours of ui 2 _ as function of stellar rotation param- 
eter D 2 and wavenumber K for a composite partial SID system 
with 8 = 0.2, /3=10 and F = 0.8. 




Figure 10. Contours of as function of stellar rotation param- 
eter D 2 and wavenumber K for a composite partial SID system 
with 8 = 0.2, /3 = 10 and F = 0.1. 



our own Galaxy, the Milky Way, for which the necessary ob- 
servational data are available, such as the gas fraction, the 
stellar and the gaseous velocity dispersions and the epicyclic 
frequency, and so forth. 

There are several inadequacies of our model. First, the 
isothermal assumption implying constant a s and a g through- 
out the disc system is a gross simplification; the velocity dis- 
persion a s in the Milky Way decreases with increasing radius 
(e.g., Lewis & Freeman 1989). Secondly, by the polytropic 
approximation, the surface mass densities are characterized 
by power-law tx r _1 profiles, which appear not to be the case 
of the Milky Way (e.g., Caldwell & Ostriker 1981). Neverthe- 
less, to estimate local stability properties of a disc portion, 
we may take a s = 50 km s _1 , V s = 220 km s _1 , <5 = 0.1 
and fiat rotation speed V s = 220 km s _1 around 4 kpc from 
the center of our Galaxy. For a massive dark matter halo 
such that F = 0.1, the relevant parameters to determine the 
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Figure 11. Contours of uP_ as function of stellar rotation param- 
eter D 2 and wavenumber K for a composite partial SID system 
with 5 = 1, /3 = 30 and F = 0.5. 



effective Q parameter for a composite partial SID system 
are S = 0.1, (3 = 50, D 2 S = 20 and F = 0.1. Following the 
definition of Jog (1996), the value of Qj > 1. While for full 
SIDs with F = 1, we get Qj = 0.17 which means local in- 
stability. Obviously, the dark matter halo plays the key role 
in stabilizing the system. 

The composite partial SID model with power-law oc r _1 
surface mass densities may not describe the Milky Way well. 
But for those disc galaxies with flat rotation curves and 
approximate power-law oc r _1 surface mass densities, our 
model should offer a sensible local criterion for axisymmetric 
instability. For example, given radial variations of velocity 
dispersions and gas mass fraction, we may deduce locally 
unstable zones. 

For circumnuclear discs around nuclei of galaxies, the 
stellar velocity dispersion a s can be as high as several hun- 
dred kilometers per second along the line of sight (e.g., Whit- 
more et al. 1985; McElroy 1994); this may then lead to small 
/)s . Within such a parameter regime, the collapse instability 
may set in during a certain phase of system evolution. For 
instance, in the case of Fig. 2 with S = 0.2 and f3 = 10, we 
choose V s = 200 km s^ 1 and a s — 300 km s _1 in a circum- 
nuclear disc region of a disc galaxy. The resulting D 2 = 0.44 
would give rise to collapse and the composite SID system 
might eventually evolve into a bulge of high velocity disper- 
sion. 

Another possible collapse situation of interest corre- 
sponds to small D 2 for sufficiently low disc rotation speed V. 
Potential applications include but not limited to protostellar 
discs in the context of star formation. We here briefly dis- 
cuss several qualitative aspects. For a protostellar disc sys- 
tem, the two disc components here may be identified with 
the relatively "hot" gas disc and the relatively "cool" dust 
disc. Interactions of radiation fields from the central proto- 
star with cool dusts in the disc lead to infrared emissions. In 
our idealized treatment, the gas and dust discs are treated 



as two gravitationally coupled thin SIDs. This simple model 
treatment may need to be complemented by other processes 
of coupling such as collisions between gas and dust particles. 
When such a composite disc system rotates at a low speed 
(e.g., Shu et al. 1987), the rotation parameter D could be 
low enough to initiate large-scale collapse. The subsequent 
dynamical processes will drive the composite SID system to 
more violent star formation activities. 

So far we have investigated the axisymmetric stability 
problem for a composite system of gravitationally coupled 
stellar and gaseous SIDs using the WKBJ approximation 
and we now summarize the results of our analysis. 

First, because the single SID case studied by Shu et al. 
(2000) may be regarded as the special case of a compos- 
ite SID system, the results of our WKBJ analysis clearly 
support the physical interpretations of Shu et al. (2000) for 
the marginal stability curves of axisymmetric perturbations. 
The "ring fragmentation regime" is related to the familiar Q 
parameter (Safronov 1960; Toomre 1964). The presence of 
the Jeans "collapse regime" can be traced to the SID model 
where self-gravity, effective pressure, and SID rotation com- 
pete with each other. 

Secondly, the recent study of Lou & Shen (2003) gener- 
alizes that of Shu et al. (2000) and describes exact pertur- 
bation solutions in a composite SID system. We obtained 
marginal stability curves for axisymmetric perturbations 
that are qualitatively similar to those of Shu et al. (2000) 
and that depend on additional dimensionless parameters. It 
is natural to extend the interpretations of Shu et al. (2000) . 
Through the WKBJ analysis here, we now firmly establish 
the presence of two regimes of "ring fragmentation" and 
"collapse" in a composite SID system. It is fairly straight- 
forward to apply our exact D— criterion for a composite SID 
system, not only for disc galaxies, but also for other disc 
systems including circumnuclear discs and protostellar discs 
etc. 

Thirdly, in the WKBJ approximation, it is also possi- 
ble to relate our D— criterion to the Qe criterion similar to 
that of Elmegreen (1995). Because k s and K g are different 
in general, our WKBJ treatment is not the same as that of 
Elmegreen (1995). 

Fourthly, also in the WKBJ approximation, we relate 
our D— criterion to the Qj criterion as defined by Jog (1996) 
but with k s / Kg in general. We find that the Qj criterion 
is fairly robust in the WKBJ regime. 

Finally, we further consider the axisymmetric stability 
of a composite partial SID system to include the gravita- 
tional effect from an axisymmetric dark matter halo. The 
stabilizing effect of the dark matter halo is apparent. 
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